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Abstract 

We propose a method to study the magnetic properties of a disordered Ising kagome lattice. The model considers 
small spin clusters with infinite-range disordered couplings and short-range ferromagnetic (FE) or antiferromagnetic 
interactions. The correlated cluster mean-field theory is used to obtain an effective single-cluster problem. A finite 
disorder intensity in FE kagome lattice introduces a cluster spin-glass (CSG) phase. Nevertheless, an infinitesimal 
disorder stabilizes the CSG behavior in the geometrically frustrated kagome system. Entropy, magnetic susceptibility 
and spin-spin correlation are used to describe the interplay between disorder and geometric frustration (GE). We find 
that GE plays an important role in the low-disorder CSG phase. However, the increase of disorder can rule out the 
effect of GE. 


1. Introduction 


Erustration in magnetism refers to competing exchange interactions that drive magnetic moments to an utterly 
conflicting situation lllyl]. There are two well known sources of frustration. The first one is caused by the geometry 
of some lattices combined with antiferromagnetic (AE) interactions in such a way that the ground state becomes highly 
degenerate without any conventional long-range order (LRO) |13l- This can give rise to the spin-liquid state, in which 
the magnetic moments are extremely correlated but still fluctuate at lower temperatures H. The second one is caused 
by the presence of disorder. Eor instance, it can produce competition between ferromagnetic (EE) and AE interactions 
arising a spin-glass (SG) state, where the magnetic moments are frozen in random directions j^]. Although geometric 
frustration (GE) is not required to stabilize the SG phase, there are several geometrically frustrated magnets showing 
SG-like behavior 111] . Intriguinly, in some of these materials the SG-like behavior appears at very low levels of disorder 
or even for what seems a disorder-free situation QlSIli . These rather unexpected results raise issues such as how 
a SG-like state can be stabilized in such conditions. 

It should be remarked that there are few theoretical works that investigate how a SG-like state can be stabilised 
in geometrically frustrated systems at very low levels of disorder SEllil fT3.fl3l] . In fact, one of the difficulties 
seems to be that the conceptual and mathematical framework to describe the disordered-driven SG state needs to be 
adapted to describe this new class of problem. Recently, an important step forward has been made by Andreanov and 
collaborators ESi. They formulated the problem in terms of a geometrically frustrated Edwards-Anderson (EA) type 
model. The disorder was introduced as a random deviation of the short-range uniform AE interactions. They showed 
that in this approach geometric frustration is quite sensitive to SG instability even for low disorder. 

Recently, a cluster spin-glass (CSG) phase 0 [II [II Hi] and short -range spin correlations within the SG-like 
phase 1 13,13 lil have been found in several geometrically frustrated magnets. Eor example, the presence of spin clus¬ 
ters can favor a SG-like phase in Zn3V308 ImI . CaBaEe 4 _jLij :07 113] and Nd 5 Ge 3 llfili . In particular, for ZnCr 204 , 
there is the suggestion that small spin clusters rather than canonical spins can bring hypersensitivity to disorder 13 ■ 
In addition, the kagome Ising antiferromagnet Co 3 Mg(OH) 6 Cl 2 provides a SG-like behavior with short-range corre¬ 
lations, apparently without disorder OlSIl . In this magnet the neutron diffraction and muon spin rotation/relaxation 
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results could also support the presence of small spin clusters. In these systems the freezing behavior may not be well 
described as a conventional SG. From the theoretical side, novel techniques are then needed to account for the cluster 
SG state found in geometrically frustrated systems. 

More recently, new approaches have been proposed using clusters of spins instead of canonical ones 


21112211 . 


Some of these theories (see, for instance. Refs. 112 iL 12211 ') follow the procedure introduced by Soukoulis and Levin 
123. 3]- The Soukoulis-Levin theory has been proposed as an improvement of the mean-field description of the 
SG problem by introducing intracluster degrees of freedom. Their Hamiltonian has two terms, the first one gives 
the intracluster short-range spin interactions which can be AF or FE. In this term, details of the original lattice can 
be included in the cluster geometry. The disorder is given by a long-range random intercluster (in fact, an infinite 
range) interaction from which a cluster SG state can appear. However, the procedure to describe the AF/FE LRO 
in the Soukoulis-Levin theory is rather artificial. The AE/EE LRO is obtained by displacing the distribution of the 
long-range random intercluster interactions by a uniform term, which is completely unrelated to intracluster AE/EE 
short-range interactions. Eor a geometrically frustrated lattice such procedure implies that GE effects coming from 
intracluster spin correlations are poorly connected with intercluster spin correlations. In view of that, it is crucial to 
account better for the role of AE interactions in this particular type of disordered lattice. 

The purpose of the present work is to study how a CSG phase can be stabilised in the disordered kagome lattice 
Ising model for low levels of disorder. The kagome lattice is one of the most straightforward examples of highly geo¬ 
metrically frustrated system being a possible realization of exotic magnetic states 3 2 ^- We follow a procedure 
similar to that of the Soukoulis-Levin theory. The lattice is divided in a set of clusters with short-range AE/EE intr¬ 
aclusters interactions. Nevertheless, the intercluster term does not have only the long-range random interactions but 
also short-range AE/EE interactions. Here, we adopt the correlated cluster mean-field (CCME) theory to treat these 
short-range interactions ll27n . In particular, the Soukoulis-Levin theory combined with the CCME method have been 
used in the disordered square lattice Ising model showing that a LRO AE evolves directly from the AE short-range 
interactions below the Neel temperature (see Ref. 1^ '). In addition, at very low temperature, these short-range AE 
interactions produce a ^n cluster compensation which makes less effective the random intercluster interaction to 
stabilize a CSG phase 12811 . On the contrary, for the disordered kagome lattice, it is expected that GE effects coming 
from the cluster core should undermine the AE LRO even at very low temperature. However, they also could en¬ 
hance the uncompensation of the spin cluster, which makes the random intercluster interaction more effective. As a 
consequence, the level of disorder which stabilizes the CSG would be lowered. In order to check how effective this 
mechanism is, we also analyse the problem when the AE interactions are replaced by EE ones. In this case, a EE LRO 
can appear, which should change the level of disorder in which the CSG phase is stabilized as compared with the AE 
case. 

The CCME theory has been initially proposed as a procedure to improve the mean-field app roximation for canon¬ 
ical spins taking into account spin correlations neglected by other mean field approaches 1271] . In the CCME theory 
the lattice is divided into symmetric clusters which are, in fact, just an artefact to incorporate spins correlations. The 
internal field acting on a given cluster is strongly dependent on the surface spin configurations of this cluster. The 
CCME theory is able to improve results for critical temperature, spin-spin correlations and critical exponents in EE 
spin systems as compared, for instance, with usual mean-field approximations 3. It is important to remark that 
in the present work we have verified whether the CCME theory applied in the clean AE kagome lattice Ising model 
can correctly capture GE effects. The results are in good agrement with the exact ones (see Ref. 1^ 1. such as absence 
of LRO at finite temperatures and residual entropy for the geometrically frutrated case. 

We also adopt the infinite-range intercluster random interaction given as the van Hemmen (vH) model (m, 321. 
The original vH model for canonical spins belongs to a class of separable on-site SG model. It is a weakly correlated 
disorder that introduces a small amount of frustrated configurations and, most importantly, it can be treated avoiding 
the replica method. This model, even lacking the multi plic ity of metastable states found in the SK model, gives 
several static quantities in good agreement with real SG 113 iL 13211 . In particular, it allows accessing thermodynamic 
quantities, such as entropy, inside the SG phase without trouble caused by the replicas. It should be mentioned that 
the vH disordered cluster model in a square lattice treated with the CCME method has been studied in Ref. 13 ■ This 
allows to compare two distinct lattice geometries, one of them with GE, both solved with the same method. 

The remainder of the paper is organized as follows. The disordered kagome lattice in a CCME approach is 
introduced in Sec. |2] Numerical results for the kagome lattice without disorder are presented in Sec. 13.11 The 
interplay between disorder and GE is discussed in Sec. 13.21 Remarks and conclusion are given in Sec. |4] 
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2. Model 


We start from the Ising model H - - 2, y JijCTiCTj on the kagome lattice, where Jij — Jq + SJij is the exchange 
interaction with a random deviation dJjj. The lattice is divided into Nd star-shaped clusters with - 12 sites (see 
Fig. ([T])), resulting in 

Vrf n, Nd 

- Z Z^-^o'' + ^-^-7 (1) 

V i.j v,A i,j 

where cr^ = +1 is the Ising spin of site i of cluster v. The first term of Eq. ([T]i represents intracluster interactions, 
in which we consider only FE or AF short-range interactions, 7^*' = Jq, without random deviation {6JJJ - 0). For 
the second one, we assume two types of couplings: the Jq between nearest-neighbor spins of neighbor clusters and 
long-range disordered couplings among all spin pairs of different clusters (d/T* = 57’'*). Explicitly, the Hamiltonian 
is represented by: 

Nd n., Nd 

y (ij) (Vi.-tj) v,A 


where (v,, Ty) denotes sum over nearest-neighbor spins of different clusters, and Sy - 2"' o-y. represents the total 
magnetic moment of cluster v. In other way, model (|2]i is a kind of hybrid between a cluster approach (CCMF) to 
the kagome lattice, which allows to take into account short-range interactions, and the introduction of disordered 
couplings between these clusters, which could induce a SG behavior. In this cluster formalism the lattice structure is 
preserved inside the clusters, in which the short-range interactions remain relevant to determine the system behavior, 
mainly in the geometrically frustrated case. In particular, the intercluster disorder follows the van Hemmen coupling 
11321] : dJ'’^ - where 7 is the strength of disorder, and ^v’s and y/^’s are independent random variables 

with Gaussian distributions of variance one. 

The intercluster disorder is evaluated as in Ref. J^. The last term of Eq. dJI is written in quadratic forms. The 
partition function for a set of quenched variables {^, j]] can then be obtained as 


^Tre 


-mm ^ J" 


dq\dq 2 dq^ exp I -N 


pj{q 




1 


- ^lnTrexp(-/3//^//) 


(3) 


where q\{^, q), qii^, q), and qj,(^, q) are introduced by Hubbard-Stratonovich transformations, N — Ndiis and Hejf - 
+ riy)qi - - ^vq 2 ]Sy - Jq I,(ij)ey (ricr j) - Y^iv^Aj) JQO-y,o-y.. In the thermodynamic limit (Nd <»), 

the functional integrals over {qni^,q)} (n - 1,2 and 3) are computed by the steepest descent method. As results 
73 = 7i + 72 and q - qi - q 2 , in which the SG order parameter q is obtained from 1132113 111 
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N \ Tr 
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where 


Nd 

Heff = ^[-7(fy -H qy)qSy- Y -^ocriCTj] - Y Joo-y^o-y. 

V {i,j)ev {vi,Aj) 


(5) 


and represents the average over ^ and q. In this approach, the disorder determines the cluster-like behavior, 

in which the cluster magnetic moment is under an effective internal field that depends on the SG order parameter q 
(see first term of Eq. (|5]l). In particular, q is an order parameter that represents the behavior of all cluster spins, which 
characterizes a CSG phase. 

In Eq. (|5]l there are still short-range couplings (FE or AF) between neighbor clusters. These spin interactions of 
different clusters are treated in the framework of the CCMF approach 12711 . which results in the following effective 
one-cluster model 


Heff ^ - Y ~ Y - Y 

(/,7)gv /gv /gv 


( 6 ) 


3 















Figure 1: Schematic representation for a kagome lattice divided into clusters with tis = 12. The mean fields are pointed by arrows that represent 
(a) the interactions between cluster v and its neighbors; (b) the interactions on v' used to evaluate the . 


where is an effective field acting on the spin cr, of site i, which belongs to the cluster border 

(v). The indices j and k represent sites of v and nearest-neighbors of site i. Hereafter, we remove the spin cluster 
index v. In the CCMF routine the mean fields and are dependent of the spin configurations of cluster v. 
For each spin configuration, {<Ti<Tj] - {++, H—, —h, —}, a different mean field is obtained: m'^~, m . For 

example, the set of mean fields } acting on spin crs is obtained from the average value of cr 2 ' of cluster v' for 

all possible configurations of the spins 0-5 and cth of cluster v (see Fig. [T]i. This average is evaluated from 


/Trcr2'e 

\ Tr.-^«W L’ 


(7) 


where s and s represent the spin states 0-5 and cri j, respectively. The Hamiltonian of cluster v' is given by 
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where the mean fields and are replaced by s and s, respectively. Here we can explore the symmetry 

of the spin cluster border to obtain all sets of mean fields. We checked the validity of this symmetry by explicitly 
evaluating all mean fields for all border spins (without any symmetry). The result was the same as that obtained by 
adopting the symmetry. 

The effective one-cluster problem is solved by computing Eqs. (HJi, ([Til and (| 8 ]l self-consistently. After, we can 
derive other thermodynamic quantities. For instance, the magnetization per site and the magnetic susceptibility are 
calculated as 


m — — 
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1 /Tr2^cr,e 
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(9) 


andx - (ff )^_Q- We also compute the staggered magnetization of three sublattice to identify a possible AF order.The 
internal energy per site is obtained from 


_ Jo / 7'r 2(/j)ev 




n.s \ 


Tre--®"'// 


7 , 




1 I TxY.ieyhf(Tie~P^‘ff 
2«s \ Tre~P^‘ff 


( 10 ) 


where the factor 1 /2 accounts for the fact that each effective interaction between clusters is shared. When 7 or ^ is 
equal to zero, the cluster is an artefact to introduce correlations and the better expression for the internal energy is 


3 rtr 
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where the factor 4/3 accounts for the intercluster interactions [331]. In particular, the 


first term between brackets of Eq. (fTOl i corresponds to the intracluster spin-spin correlation function {cricrj)intra- The 
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Figure 2: Entropy and magnetic susceptibility as a function of the temperature for the geometrically frustrated kagome lattice without disorder, in 
which s(T —> 0) = 0.503. The inset exhibits the behavior of magnetization, and s for ferromagnetic interactions. 


specific heat (Cy) is obtained by numerical derivation of the internal energy. The entropy per site is given by integration 
of Cv/T: 


s — 





( 11 ) 


3. Results 

Numerical results are obtained from the exact diagonalization of the effective single-cluster problem. We obtain 
other thermodynamic quantities such as magnetization, magnetic susceptibility, entropy, and spin-spin correlation 
function. We also build phase diagrams of the temperature T versus disorder intensity J. The results are in unity of 
l/oL which for numerical purposes is used I/qI = 1- The cluster spin-glass (CSG) phase occurs when q and m - 0, 
and the ferromagnetic (FE) order is characterized by m + 0 and q -Q. The antiferromagnetic order is not found. 


3.1. Kagome lattice within the CCMF 

Figure [^presents the results for the pure kagome lattice {J - 0). It exhibits the entropy s and magnetic suscep¬ 
tibility for the geometrically frustrated case (7o < 0). The entropy curve goes towards a non-zero value as T ^0: 
S(T —> 0) = 0.503. This residual entropy is a consequence of the strong GF and indicates a high ground-state degen¬ 
eracy. This finding is in perfect agreement with the exact result for the AF Ising kagome lattice ll34ll . Furthermore, 
the magnetic susceptibility diverges as T —» 0; however, there is no phase transition at finite temperature S. The 
residual entropy and the absence of FRO suggest a spin-liquid ground state |01 that is in agreement with early results 
1 3511 . Therefore, the CCMF theory arises as a good alternative to study this strongly geometrically frustrated system. 
In addition, the inset of Fig. |2]shows results for thjae FE kagome lattice (/q > 0). The magnetization and susceptibil¬ 
ity indicate a para mag netic PM-FE phase transition at Tc - 2.417o that overestimates around 12% the exact solution 
(j-exact _ 2.147o 13411 '). The entropy also exhibits a typical curve for the Ising FE kagome lattice. It is important 
to remark that the CCMF theory applied to the kagome lattice provides better results (as compared with the exact 
solution) than conventional mean-field theories 


3.2. Kagome lattice with disordered clusters 

The effect of disorder on the magnetic properties of the kagome lattice is now analyzed. For instance, the FE 
long-range order is found in the case without GF (Jq > 0). Figure ©a) shows that the FE order is still stable for low 
intensities of disorder. However, the CSG behavior appears when J is higher than a critical value. There is a range 
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Figure 3: Phase diagrams of the critical temperatures as a function of the disorder intensity, (a) FeiTomagnetic interactions, (b) Antiferromagnetic 
interactions. The inset exhibits the Tf behavior for high disorder intensities. 


of disorder in which a reentrant PM-CSG-FE phase transition occurs when the temperature decreases. In particular, 
the CSG-FE transition is first-order and it is located by comparing free energies obtained from: F - U - TS. To 
summarize the CSG phase occurs only at high disordered regimes in the kagome lattice without GE. 

Figure (13-b) exhibits Tf as a function of J for the AF kagome lattice, in which the interplay between GF and 
disorder brings on a completely different scenario. There is no FRO at finite temperature for 7 = 0, but an infinitesimal 
disorder is enough to change the spin-liquid ground state into the CSG phase. It means that the GF helps introducing 
the CSG phase at very small intensities of disorder. Furthermore, there are clearly two distinct behaviors for Tf{J). 
One for weak disorder, in which the effect of GF is relevant. Other for high intensities of disorder, where the Tf(J) 
inclination is more pronounced (see inset of Fig. [3-b). 

In the following, the interplay between GF and disorder is discussed in terms of some thermodynamic quantities. 
For instance, the entropic behavior of the AF kagome lattice is strongly affected by disorder. Figure |4] shows that 
a very low disorder intensity leads the entropy to zero when T ^ 0. This indicates that an infinitesimal disorder 
is able to vanish the residual entropy and stabilize the CSG phase. It is important to point that the CSG phase is 
between clusters. In fact, the system ground state becomes cluster CSG. In other words, the entropy is an extensive 
quantity that characterizes the whole system, but it can not identify spin states inside the clusters. For higher disorder 
intensities, the entropy also decreases towards zero when the temperature diminishes (see the inset plot of Fig. |4]i. 

Figure|5]-a shows the magnetic susceptibility, which presents a flat plateau below Tf that is a characteristic of spin 
glass phase. For high temperatures, the system follows the Curie-Weiss law. These features are found for any finite 
disorder. The linear extrapolation of;^f * (see the inset in Fig. |5}a) allows to estimate the Curie-Weiss temperature 
Ocw ~ -5|7o| that is independent of 7. In the present approach, the;^^ of the PM phase is not affected by the disorder 
intensity (see Fig. |5}a). As a consequence 9cw becomes independent of 7. However, an important result is obtained 
from the frustration parameter fp — \6cw\IT / when the disorder changes. For a weak disorder regime, the frustration 
parameter is higher than for the strong one: e.g., for 7/|7o| = 0.5 {Tj - 0.167|7o|), fp ~ 30, and for 7/|7o| - fp ~ 3. 
It means that the GF remains extremely important for weak disorders, but its influence is quite reduced for stronger 
disorders. 

Other result that indicates the presence of two disordered regimes is the correlation function (cr,cry),„,ra (Fig. 
13-b). This intracluster correlation function considers only interactions between spins of the same cluster. It gives 
information on the spin-spin correlation inside the clusters, but not among different clusters or on the whole system. 
Figure 13-b shows that this correlation is negative for weak disorders. In particular, {cricrj)j„,ra is very close to that 
one for the pure AF kagome (see inset of Fig. |5}b). However, it becomes positive at low temperatures for stronger 
disorders. This change in the correlation function suggests that pairs coupled ferromagnetic ally inside the cluster are 
energetically favored by disorder. As a consequence, the cluster magnetic moment increases, which can favor the 
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Figure 4: Entropy as a function of the temperature for the increase of disorder intensity in the AF case. The inset shows a result for a high value of 
J. 


intercluster disordered couplings. These results indicate that a high intercluster disorder overcomes the intracluster 
interactions, supressing the effects of GF. 




Figure 5: (a) Magnetic susceptibility versus temperature for several values of J. The inset exhibits the inverse of susceptibility and the dotted line 
stands for the linear extrapolation from high temperatures, (b) Spin-spin coiTelation function behavior for different intensities of disorder. The inset 
shows results for lower J values. 


4. Remarks and Conclusion 

We study the interplay between GF and disorder by using the Ising kagome lattice in a cluster method. The 
model adopts a van Hemmen type of disordered interaction between clusters combined with the short-range CCMF 
approach. This treatment allows us to access thermodynamic information of the disordered kagome lattice and draw 
phase diagrams to analyze the competition between CSG and conventional orders. 

The CCMF method improves the mean-field description of the pure kagome lattice. In particular, the critical 
temperature for the FE kagome lattice and the absence of long-range order in the AF case are in good agreement 
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with the exact results 34]. The thermodynamic description of the entropy, magnetic susceptibility and spin-spin 
correlation function indicates that the CCMF theory is a valuable tool to study geometrically fmstrated systems. 

The presence of disorder in the FE kagome lattice introduces competition between ferromagnetism and CSG 
state. The FE order is found at low intensities of disorder. However, the CSG phase appears at sufficiently strong 
disorder. The disorder actives the presence of clusters, in which the intracluster EE interactions introduce a high 
cluster magnetic moment. This situation enhances the effect of the intercluster disorder that favors the CSG behavior. 
It means that a critical value of disorder is needed to establish the CSG phase. 

On the contrary, an infinitesimal disorder is able to stabilize the CSG behavior in the geometrically frustrated 
kagome lattice. In this case the lattice is divided into clusters for any disorder intensity. The entropy goes to zero 
when the temperature decreases and there is no residual entropy. It suggests that any amount of disorder is enough to 
stabilize the CSG phase instead of the spin-liquid ground state. Nevertheless, the frustration parameter is extremely 
high for weak disorders. These results indicate that the GE of the whole lattice is replaced by geometrically frustrated 
clusters, which helps to introduce the CSG ground state for an infinitesimal disorder. 

This interpretation is also supported by the freezing temperature as a function of disorder strenght that shows 
two different behaviors. Eor weak disorders, the Tf{J) inclination is smaller than in a strong disorder regime. The 
intracluster spin correlation function indicates that the AE interactions remain important inside the cluster for weak 
disorders. These frustrated intracluster interactions favor an uncompensated cluster magnetic moment, which is small 
but sufficient to stabilize the CSG behavior at weak intensities of J. As disorder increases, the intracluster AE inter¬ 
actions lose importance and the spin-spin correlation becomes positive. Thus, a higher cluster magnetic moment is 
energetically favored by the intercluster disordered interactions. It increases the Tf(J) inclination at the same time that 
the frustration parameter becomes very small. In other words, an infinitesimal disorder is enough to bring a cluster 
CSG phase in which the GE plays an important role. However, as long as disorder increases, the effect of GE is ruled 
out. Our results suggest that the formation of very small spin clusters can contribute to stabilize a spin glass behavior 
even at an extremely low level of intercluster disorder. In this context, the GE helps to introduce uncompensated 
spin clusters, which become hypersensitive to disorder. Although we study a particular model, the results for weak 
disorders can shed some light in the unconventional SG-like behavior in the compound Co 3 Mg(OH) 6 Cl 2 ra . In 
particular, the results reported in figures 4 and 5 (a) for entropy and magnetic susceptibility are qualitatively similar 
to those found in this compound. 
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